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Abstract 

Building on an analogy between the ageing behaviour of magnetic systems and grow¬ 
ing interfaces, the Arcetri model, a new exactly solvable model for growing interfaces is 
introduced, which shares many properties with the kinetic spherical model. The long-time 
behaviour of the interface width and of the two-time correlators and responses is analysed. 
For all dimensions d ^ 2, universal characteristics distinguish the Arcetri model from the 
Edwards-Wilkinson model, although for d > 2 all stationary and non-equilibrium expo¬ 
nents are the same. For d = 1 dimensions, the Arcetri model is equivalent to the p = 2 
spherical spin glass. For 2 < d < 4 dimensions, its relaxation properties are related to the 
ones of a particle-reaction model, namely a bosonic variant of the diffusive pair-contact 
process. The global persistence exponent is also derived. 
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1 Introduction 


The physics of the growth of interfaces is a paradigmatic example of the emergence of non- 
eqnilibrium cooperative phenomena [T, ,21 El 13 Gb HE I]E The morphological evolution of the 
growing interface is characterised by its fractal properties as well as several growth and rough¬ 
ness exponents. More recently, aspects of the time-dependent evolution of the interface have 
been investigated and have been found to be quite analogous to phenomena encountered in 
the physical ageing in glassy and non-glassy systems mm- It has been understood that grow¬ 
ing interfaces can be cast into distinct universality classes, including those associated to the 
Edwards-Wilkinson (ew) equation [TO] , where the relaxation is dominated by the linear in¬ 
terface tension, and the Kardar-Parisi-Zhang (kpz) equation m, which includes as well a 
non-linear effect of the local slope. 

The exponents of growing interfaces are defined as follows. The basic quantity is the time- 
space-dependent local height h(t, r ) and one is mainly interested in the fluctuations around the 
spatially averaged height h[t) := L~ d J2 re \h(t,r)] for notational convenience defined here on 
a hyper-cubic lattice A C with L d sites and the sum runs over the lattice sites. A central 
quantity is the interface width 

w 1 2 (t;L) := jdJ2({h(t,r)-h(t))) 2 = L 2a f w (tL~ z ) ~ { ^ 2q ; J \ (1-1) 

and we recall the generically expected Family-Viscek scaling form [12], where a is the rough¬ 
ness exponent, /3 the growth exponent and 2 = a/(3 the dynamical exponent. Herein, (.) 
denotes an average over many independent samples (under the same thermodynamic condi¬ 
tions). Throughout this work, we shall consider the L —> oo limit. The initial state will always 
be a flat, uncorrelated substrate. Relaxational properties of the interface can be characterised 
by the two-time correlations and (linear) responses 

) (h(s, 0) - (h(s)))) = s- b F c (1.2) 

= (h{t, r)h(s, 0)^ = s~^ a F R 
j =o 

where spatial translation-invariance has been implicitly admitted and j is an external field 
conjugate to h. Furthermore, in the context of Janssen-de Dominicis theory, h is the conjugate 
response field to h USE]- In the long-time limit, where both t,s r micro and i-s > Tmicro 
('hnicro is a microscopic reference time), generalised Family-Viscek scaling forms were assumed 
mmmmm which are analogous to the scaling forms of ‘simple ageing’ in other non- 
equilibrium systems [8[IT9ll9]. The exponent b = —2/3 [13, ED], but the relationship of a to other 
exponents seems to depend on the universality class pS]- Finally, the autocorrelation exponent 
Ac and the autoresponse exponent X R are defined from the asymptotics Fc,r(i/ : 0) ~ y~ x c,RF 
as y —> oo. There is a bound Ac > (d + zb)/ 2, see appendix B. Field-theoretical considerations 
in ageing simple magnets with a non-conserved model-A dynamics and with disordered initial 
conditions strongly indicate that the non-equilibrium exponents Ac, X R should be independent 

1 Applications include settings as distinct as deposition of atoms on a surface, solidification, flame propagation, 

population dynamics, crack propagation, chemical reaction fronts or the growth of cell colonies. 



C(t,s;r) ■= ((h(t,r)-(h(t)) 
6 (h(t, r) — h(t)) 


R(t,s\r) := 


Sj(s, 0) 
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of those describing the stationary state [T3] 2T| [7j. This is different, however, for the KPZ 
universality class, where for dimensions d < 2 it was shown that Ac = d , to all orders in 
perturbation-theory |22j . However, for d > 2 there exists a strong-coupling fixed point of the 
KPZ equation which cannot be reached by a perturbative analysis, see e.g. [23j|24j[7]; but which 
can be studied through non-perturbative renormalisation-group techniques [25l [26] . 

Known estimates of all these exponents, for several universality classes, are listed in table [Q 
In most of the quoted experiments, two-time autocorrelators such as C(f, 0; 0) were also used 
to extract (3, these approaches rely on the relation b = — 2/3@ The study of two-time global 
quantities, also of interest for experiments, gives a different access to the exponents, especially 
/3 [4TJ 12] • For the experiments reported in [36] [33] [32], merely an average of the values 
quoted in the sources is included in table [Q While the theoretical analysis, at least for the 
exponents a, /?, z in the EW and KPZ classes, has been achieved long ago mm, recent advances 
in experimental techniques have permitted to obtain precise and reliable estimates of these 
exponents, independently and in several distinct systemsjl and in good agreement with each 
other and with the theoretical predictions^ Notably, the great variety of systems for which the 
exponents z, /3, a are consistent with the KPZ-equation, clearly attest its universality, at least for 
ID systems. For the ID (positive) quenched KPZ class, a mapping onto ID directed percolation 
was proposed [3Ej- This produces a = z/j_/zzy — 0.63261 which apparently fits the experimental 
and simulational data well (using the very precisely known values of the exponents of 
directed percolation, see e.g. BSD- However, other predictions following from this mapping, 
such as z = 1 [38] [I], do not seem to be generically reproduced by the presently available 
evidence. In several of these experiments, universal amplitude ratios have been measured as 
well [371 US] [32] |35]. For further details on procedures and the extraction of the exponents, we 
refer to the quoted sources. For a recent thorough review on experimental results, see m- 

In d — 1 dimension, several further notable advances have been achieved. These concern a 
remarkable exact solution of the KPZ equation and the spectacular relationship of the proba¬ 
bility distribution V{h) of the fluctuation h—h with the extremal value statistics of the largest 
eigenvalue of random matrices, see m m edi eh eh ng, and its successful experimental con¬ 
firmation da m 03 m izhi gsj. Still, this progress seems to rely on specific properties of 
the one-dimensional case. One may ask, if there might exist alternative routes to a further 
conceptual understanding, less dependent on the particular mathematical circumstances which 
render the ID KPZ equation analytically treatable. Here, an analogy with magnetic systems 
at their critical point might be helpful. Some elements of such an analogy are sketched in ta¬ 
ble [2] Indeed, several observables derived from the main quantity, namely the order parameter 
4>(t,r ) or the interface height h(t,r), respectively, show an analogous dynamical behaviour. In 
particular, the exponent b, to be read off from the scaling behaviour (ll.2jl of the autocorrela¬ 
tor, is b = 2 j3/{vz) for a critical magnet and b = —2/3 for interfaces. This is consistent with 
the scaling behaviour of the equal-time variance/interface width, respectively. On the other 
hand, for magnets relaxing towards an equilibrium state, a fluctuation-dissipation theorem gives 
a = b = 2 /3/(iyz), see e.g. [9], whereas no such generic result is known for the fundamentally 
non-equilibrium interfaces ]j The analogy between magnetic systems and interfaces is further 

2 However, this relation needs no longer to hold true when d > d*. 

3 In [33], colonies of both benign and malign cells, for both flat and circular interfaces were measured. 

4 See 03]HUES] for recent estimates, of the exponents a, /3 and also of universal amplitude ratios, in the KPZ 
universality class, up to d = 11. 

5 In the EW model, one has a = b mg, while in the ID KPZ model l + a = b + 2/z 0] and the relationship 
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Tabic 1: Exponents of growing interfaces in four universality classes, namely Kardar-Parisi-Zhang (kpz), (positive) quenched KPZ 
(qkpz), Edwards-Wilkinson (ew) and Arcetri (for both T = T c and T < T c ) are represented. In each group, labelled by the universality 
class, first are indicated available exact results^ and/or recent simulational estimates.^ Then follow experimental results, where the 
system’s dimension d and the nature of the interface are indicated. The numbers in bracket give the error in the last digit(s). 

*The scaling relation a = z/3, resp. a + z = 2 for the KPZ class, was used to complete the entries as much as possible, with an error estimated from the 
relative errors of the exponent(s) given in the source. 


model 

d 

z 

0 

a 

a 

b 

Ac 


Ref. 

KPZ 

1 

3/2 

1/3 

1/2 

-1/3 

-2/3 

1 

1 

nnEan] 


2 

1.61(2)* 

0.2415(15) 

0.393(4) 

0.30(1) 

-0.483(3) 

1.97(3) 

2.04(3) 

m 


2 

1.61(2)* 

0.241(1) 

0.393(3) 


-0.483 

1.91(6) 


m 


2 

1.61(5) 

0.244(2) 

0.369(8) 





m 


2 

1.627(4)* 

0.229(6)* 

0.373(3) 





<m 

Ag electrodepositknd 

1 


** 1/3 

« 1/2 





ra 

slow paper combustion^ 

1 

1.44(12)* 

0.32(4) 

0.49(4) 





ED 

liquid crystal^ 

1 

1.34(14)* 

0.32(2) 

0.43(6) 


« -2/3 

zz 1 


ED 

liquid crystal^ 

1 

1.44(10)* 

0.334(3) 

0.48(5) 


« -2/3 

0? 


m 

cell colony growth^ 

1 

1.56(10) 

0.32(4) 

0.50(5) 





Isa m 

(almost) isotrope colloids^ 

1 


0.37(4) 

0.51(5) 





ED 

autocatalytic reaction front f 

1 

1.45(11)* 

0.34(4) 

0.50(4) 





E5) 

CdTe/Si(100) film 1- 

2 

1.61(5) 

0.24(4) 

0.39(5)* 





ED 

QKPZ 

1 


« 0.9 

0.63261 





(38113^ 

cell colony growth (disordered) ^ 

1 

0.84(5) 

0.75(5) 

0.63(4) 





ED 

autocatalytic reaction front f 

1 


0.61(5) 

0.66(4) 





ESI 

strongly ellipsoid colloids^ 

1 


0.68(5) 

0.61(2) 





ED 

EW 

< 2 

2 

(2 - d)/ 4 

(2 — d)/2 

d/2 — 1 

d/2 - 1 

d 

d 



2 

2 

0(log)# 

0(log)# 

0 

0 

2 

2 

mm 


> 2 

2 

0 

0 

d/2-1 

d/2 - 1 

d 

d 


sedimentation/electrodispersion 

2 


0(log) # 

0(log) # 





na 

Arcetri T = T c 

< 2 

2 

(2 - d)/4 

(2 — d)/2 

d/2-1 

d/2 - 1 

3d/2 - 1 

3d/2 - 1 



2 

2 

0(log) # 


0 

0 

2 

2 



> 2 

2 

0 

0 

d/2-1 

d/2 - 1 

d 

d 


T <T C 

d 

2 

1/2 

1 

d/2-1 

-1 

d/2-1 

d/2-1 



fflat interface, ^circular interface. #For d= 2, one has w(t;L ) ~ y/lnt f w (lnL/lnt). 

























Table 2: Analogies between the critical dynamics in magnets and growing interfaces. Several 
observables derived from the central physical quantity, the order parameter/height, respectively, 
are shown, where h or j, respectively, are conjugate to the order parameter/height. The 
average (.) c denotes a connected correlator. Some models, with the equilibrium hamiltonian 
for magnets, are defined through their kinetic equations, where 77 is a standard white noise, A 
the spatial laplacian and D,g,v, p are constants. The existence of known exact solutions in the 
Ising and kpz models is indicated. 


magnets 

interfaces 

order parameter 0 (f, r) 

variance ((</>(t,r) — ( 0 (f, r))) 2 ) ~ £~ 2 / 3 /OA 

autocorrelator C(t,s ) = (0(t, r)</>(s, r)) c 

autoresponse R(t,s)= 5{(p{t,r)) /5h(s,r)\ h 0 

height h(t, r ) 

width w 2 {t) = ((h(t,r) — h{t )) 2 ) ~ t 2> 3 

autocorrelator C(t,s ) = (h(t, r)h(s, r)) c 
autoresponse R(t,s ) = 5(h(t,r))/5j(s,r)\- 0 

Models: 

gaussian held "H[0] = — ^ Jdr (V</>) 2 

d t 4> = DAcj) + T] 

Ising model %[</>] = — \ Jdr [( V</>) 2 + §</> 4 ] 

d t (p = D(A</> + gcft 3 ) + 77 

EW 

d t h = uAh + 77 

KPZ 

dth = vAh + f (Vh) 2 + 77 

Exact solutions for the Ising and kpz universality 

equilibrium d = 2 [511 EH [56] 

relaxation d = 1 [57] 

1 classes: 

relaxation d — 1 [HI HHl [51] 


illustrated by the form of the kinetic equations of the time-dependent order parameter, schemat¬ 
ically written as d t 4> = —DdR/Scj) + 77 , where 77 is a gaussian white noise. Clearly, the gaussian 
held/EW model share mean-held characteristics, whereas the behaviour of the Ising and KPZ 
models is determined by the non-linearities in their equations of motion. These non-linearities 
make exact solutions so difficult to find: for instance, in the 2D Ising model at equilibrium, 
only the cases (i) of a vanishing magnetic held h = 0 and arbitrary temperature T f5T] and (ii) 
of the fixed temperature T = T c and an arbitrary magnetic held h [55, 56] have been solved 
(see [58] for a recent overview of numerical and experimental tests). 

In magnetic phase transitions, in view of these technical difficulties, the so-called spherical 
model was originally proposed by Berlin and Kac [59] in order to be able to explore generic 
properties of equilibrium phase transitions in the context of a non-trivial exactly solvable model. 
A simple way to introduce it is to replace the discrete Ising spins cq = ±1, attached to each site 
i of the lattice A, by continuous spin variables (jj 1 —> e R subject to the ‘spherical constraint’ 

A s 2 == A/", where A f is the number of sites of the lattice. Especially in the form of the 
‘mean spherical model’ (where the spherical constraint is only required on average [03], which 
considerably shortens the calculations, without modifying the critical behaviour), the spherical 
model has become an often-used test case in many distinct situations. For the most common 
case of short-ranged interactions, a phase transition with a critical temperature T c > 0 exists in 
dimensions d > 2. For dimensions 2 < d < 4, the equilibrium exponents are different from the 
mean-held theory of the free gaussian held [59], For the dynamics, and in the continuum limit, 
this amounts in the kinetic equation of motion to replace the non-linearity by a more simple 

for d > 1 is unresolved E3, see table [T] 
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term, viz. 0 3 i—>■ (0 2 ) 0 =: 3 (f)0 and to fix the Lagrange multiplier 3 (t) through the spherical 
constraint. This procedure can be applied to study the relaxational dynamics and to extract the 
exponents of both equilibrium and non-equilibrium critical dynamics, and was performed many 
times, see e.g. ( 6 H1621 |63l EH 1651 | 66 l 1671 EH [691 Eol IZH1721 ESI EH [75l |76l [771 ESI |79l |80l EH IH 2 ] 
and references therein. For reviews, see Isa ei m, uni in ej. 

Here, we shall inquire whether a spherical model variant can be sensibly defined for models 
describing growing interfaces. I 11 particular, we shall discuss the following questions: 

1 . which property of the interface heights could be construed to take only values ±1 in order 
to identify a sensible ‘spherical approximation’ ? 

2. does one obtain a non-trivial model, distinct from the mean-field-like EW class ? And 
does there exist@ a well-defined upper critical dimension d* ? 

3. in spite of the analogy with magnetic systems relaxing towards equilibrium, is the relax¬ 
ation process an equilibrium or a non-equilibrium one ? 

This work is organised as follows: section 2 defines the Arcetri model and section 3 outlines 
the exact calculation of the interface width and of the two-time responses and correlators. The 
ID model is shown to be equivalent to the p = 2 spherical spin glass. In section 4, the long-time 
behaviour is explicitly found. We discuss in detail the non-equivalence with the EW-model for 
all dimensions d 7 ^ 2 ; the relationship of the relaxation behaviour with the one of a bosonic 
particle-reaction model for dimensions 2 < d < 4; and comment on the global persistence 
probability. We conclude in section 5. Details of the calculations are treated in appendix A. 
The Yeung-Rao-Desai inequality is revisited in appendix B. 


2 The Arcetri model 

We begin by stating the definition of the modefl to be analysed, first in d — 1 dimensions. 
Consider a set of height variables h n [t) attached to the sites n of a ring with N sites. In what 
follows, we shall work with its (discrete, symmetrised) derivatives 

u n (t ) := ^ (h n+1 (t) - h n -i(t)) (2.1) 

Furthermore, to each lattice site one attaches a gaussian random variable rj n (t), with the mo¬ 
ments 

(Vn(t)) = 0 , (r} n {t)r} m ^)) = 2TT5(t - t')8 n , m (2.2) 

The defining equations of motion for the ID Arcetri model are 

d t u n {t) = V ( U n+1 (t ) + Un-iit) - 2u n (t)) + 3 (t)u n (t) + i (Vn+l (t) ~ Vn-l(t)) (2.3) 

JV-l 

^( M „(«) 2 ) = N (2,4) 

n =0 

6 The existence of a finite d* would be analogous to critical magnets, whereas for the KPZ equation the 
question is still unresolved, although numerical results suggest that d* might be infinite, see [43 ] l44l [29l . 

'The name is inspired by the lieu where this model was conceived and this work was done. 
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Figure 1: Schematic illustration of a It) growing interface respecting the RSOS constraint 
that the local slopes u i+ y 2 = h i+ 1 — h* = ±1 on nearest-neighbour sites. The interface grows 
through sequential absorption of square-shaped particles. Those ‘active’ positions into which an 
adsorption is possible are underlined in green and those ‘inactive’ ones where this is prohibited 
by the RSOS constraint are underlined in red. Active and inactive positions have to be updated 
after each adsorption process. 

where the Lagrange multiplied 3 (t) is determined from the ‘spherical constraint’ (j2.4[) and T 
and T are constants. 

Our motivation to study this model comes from the well-known representation of lattice 
realisations of the KPZ universality class in terms of the totally asymmetric exclusion model 
(tasep) HH m T\- In figure [1] an interface is sketched, and its growth occurs through the 
adsorption of additional blocks, one at a given time, but such that between neighbouring sites 
the RSOS-condition hi + \{t) — hi(t ) = ±1 is obeyed. Below the height configuration, the value of 
the interface u i+ 1/2 = h t+ i — hi on the dual lattice (the links) is indicated. Let u i+ 1/ 2 = — 1 1 —> • 
represent a dual site occupied by a particle and tq+i /2 = +1 •—> o an empty dual site. Then the 
adsorption process corresponds to an irreversible change »o —y o» of a particle hopping to the 
right between two neighbouring sites of the dual lattice, whereas the inverse reaction o« —y *o 
is not admitted. In the continuum limit, this process is described by the KPZ equation, as given 
in tabled 

In our search for a spherical model analogue for the KPZ equation, we shall therefore identify 
the slopes (12. ip as the variables on which one could impose a spherical condition (12.4p . Then, 
it should be preferable to consider the equation of motion of the slope, which in the continuum 
limit becomes u(t,x ) = d x h(t,x ) and obeys the noisy Burgers equation 

d t u(t, x) = udlu(t, x) + pu(t, x)d x u(t, x) + d x r](t , x) (2.5) 

Finally, we replace the non-linearity ud x u 1 —> 3 (t)u in (12.5p by a spherical model construction, 
analogously to the definition of the spherical model of a ferromagnet [59; [60] and fix the La¬ 
grange multiplier 3 (t) from the spherical constraint ( 1 2.4 j) Eqs. (12. 31b (12. 4p give the discretised 

8 We do not consider any fluctuations in j(t), since they do not contribute to the spatially local averages [723 
on which we concentrate here. 

9 See [86] for a rigorous proof in ID. Another often-used lattice representation is the Kim-Kosterlitz model 
[87] . where hi + 1 — hi — ±1, 0, such that the continuum limit is again the KPZ equation. 

10 In the kinetic spherical model of a ferromagnet, using a mean spherical constraint simplifies considerably 
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version of this construction. 


Clearly, eqs. fl2.4ll . fl2.5D are identical to those of the conventional spherical model, see e.g. 
[55 ., BE], but with the standard white noise rj(t, x ) replaced by its derivative d x rj(t, x). This new 
feature brings down the upper critical dimension from d* = 4 to d* = 2 [Tj. It is more surprising 
to find that the mere modification of the noise is enough to render our model identical to the 
spherical spin glass [53] . as we shall show in section 2.4. On the other hand, the present work 
with its specific way of introducing the spherical constraint is merely meant as a first step in 
a more systematic exploration of ‘spherical variants’ of equations such as (12.51) . We shall come 
back in section 5 to other possibilities which might keep a more clear memory of the non-linear 
terms in (12.51) . Higher-order non-linearities will either turn out to be irrelevant or else may 
create some ‘multicritical’ variant of scale-invariant growth EH El Kl¬ 
in order to see how to generalise this model to any dimension d ^ 1, we briefly consider the 
2D case, for simplicity of notation. The height on each site is Then we have slopes in 

each of the two lattice directions, namely 


^n,m(0 • 2 (^bi+l,ra(^) tln —1 7 ^ (^) hn,m—l (^)) (2.6) 

which obey the equations of motion (for brevity, the dependence on t is suppressed) 


dtVn.m, T (Un+l,m T U n —l t m T U n m ^i T Un,m—1 n ,m ) T <3 (^)^n,m T ^ ij]n+l,m Vn—l,m) 


dt^n,m T (Un+l,m T ^n—l,m T ^n,m+l T 1 4u n rn ) T T ( J]n,m+1 'dn,m— l) 


(2.7) 


Considering the spherical constraint, we must count how many links exist on a torus § 1 <8)S 1 with 
N x N = N 2 sites. Since on each line of the lattice, one has N links, there are N such lines and 
there two dimensions along with the links must be counted, there is a total of N x N x 2 = 2N 2 
links. Hence, the spherical constraint must be 


N-l 

£ + =2iV 2 (2.8) 

n,m =0 


The extension to generic d is now obvious and will be written down explicitly in the next 
section {2] 

Clearly, in analogy with the EW and KPZ equations, we have already performed implicitly 
a Galilei transformation and shall always work in a reference frame which moves upwards with 
the average rate of particle deposition. The parameter T in the noise correlator serves to 
distinguish the kinetic noise from the effective diffusion constant v. 

the calculations, [65] . Provided the infinite-system limit N —> 00 is taken before the large-time limit, the 
consideration of a non-averaged spherical constraint leads to the same results for the long-time behaviour [68] • 
11 The formulation chosen here is explicitly invariant under spatial rotations. Any attempt of a generalisation 
towards more than one spherical constraint will likely lead to a breaking of that symmetry. 
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3 Solution 


3.1 Equations of motion 


The solution of the Arcetri model equations of motion, i.e. eqs. (l2.7jUI2.8jl for d — 2, is 
carried out in Fourier space. In d spatial dimensions, there are d slopes, denoted u atn (t), with 
a = 1,..., d and n = (rq,..., n d ). On a torus T d = S 1 <g) • • • <g) S 1 = (S 1 ) 0 ^, represented by a 
hyper-cubic lattice with N d sites, one has the Fourier transforms 

N-l N-l 

u a (t,p ) = ^■■•^exp 

ni=0 rcd=0 

N—i N— i / 27ri \ 

u a , n (t) = ^ ^ ^exp fnj w a (t,p) (3.1) 

Pi=0 Pd=0 ' ' 


27ri , , , 

~W p n ' Ua ^> 


Our equations of motion take the form 


d t u a (t,p ) 


-2Tu(p)u a (t,p) +$(t)u a (t,p) + i sin 



Pa rf(t,P ) 


(3.2) 


with the dispersion relation a j(p) = i [l — cos {j^Pj)]’ an d ^ ie gaussian noise correlators 
(v(t,P)) = 0, {rf(t,p)rf(t',q)) = 2rTN d 5(t-t')5 p+qt0 (3.3) 

Since the definition (12. 6 p implies integrability conditions^! the interface height is given by 

u a {t,p) = isin (^rp a ^jh(t,p) ; if p ± 0, a = 1,. .., d (3.4) 

We restrict attention to interfaces which on average are initially flat and uncorrelated. Indeed, 
the initial modes h( 0 , p) are assumed gaussian random variables with the moments 

(M0,p)) = N d H 0 dp o , (h(0,p)h(0, g)) = N d H\ h p+9 , 0 (3.5) 


where H 0 i are initial parameters. The first condition (13.5p means that the average initial slopes 
vanish. From the second, if H\ — H^, the initial interface width vanishes. We shall restrict to 
this case below, see eq. (13.221) . 

The equations of motion (13. 2j) have the following solution, with a = 1 ,... ,d 


u a (t,p ) = w a ( 0 , p) exp 


-2Yu(p)t+ / dt' 3 (f) 


rt / 2 tt 


+ / dr i sin ( —p a ) fy(r, p) exp 


-2Yu(p){t - r) + / df 3 (f) 


(3.6) 


12 From (12.61) . one has in 2 D that ^(u n ,m+i{t) — u n , m -i(t)) = —v n -i im (t)), which in the continuum 

limit amounts to du/dy = dv/dx. Standard calculus [88] p. 104] then shows that there exists a function h such 
that u = dh/dx and v = dh/dy, on a simply connected domain. Eq. (13.41) expresses this in discrete Fourier 
space. 














In order to work out the explicit form of the spherical constraint, we follow the lines of the 
magnetic spherical model and define [65] 


(3.7) 


g(t) — exp f-2 J df' j(f)i 


Then the spherical constraint becomes, using (13. ip and the initial condition (I3.4p . (l3.5p . 

d TV-1 TV-1 

dN“ = E E" • E (“«.»w 2 > 

a= 1 ni=0 Tid = 0 
N—l N—l r 


E-E 

Pi=0 Pd=0 


H\X(p)e 


4Tu( P )t_}_ _|_ f dr 2 rTA(p)e^ 4raj(p)(t " T) ^ 






(3.8) 


with A (p) \= Y^ d a =\ sin 2 (jj-Pa) ■ If we now define 


1 


TV-1 TV-1 


/(*) : =NdJ2---J2 X ^' 


,- 4 Ttu(p) 


(3,9) 


P1=0 Pd=0 

the spherical constraint can be written as a Volterra, integral equation 

Hi fit) + 2TT f dr g(r)f(t - t) = dg{t) 


(3.10) 


Its solution will be found in complete analogy with well-established techniques 

At this point, the infinite-size limit N —y oo can be taken. The solution (13. 6 p becomes 


«„(t, p) = «.(0, p)e- 2r ‘"<>>> -4= + J‘dr i sin Pa p(r. MP 


(3,11) 


where the dispersion relation now reads uj(p) = X)a=i [1 — cos Pa] and p E B = [—7r, 7r] d is in 
the Brillouin zone. Letting now A (p) := sin 2 p a , we have from (3.9), 

1 


m = 


dp X(p) e 


—AYtu){jp) 


= d 


(2n) d 

— [ 

27 Tj_ 

e~ AVt Ii(AYt) . _ 


dp sin 2 p e 4Ttcosp 


4 Ft 


(e- 4 r %(4n)) 


h.L 

d -1 


dg e 


d —1 

4rtcosg ^ 


(3.12) 


where the J n denote modified Bessel functions [89]. The function g(t) is found from the Volterra 
integral equation (13. 10p . Furthermore, since (13. 4p is now written as u a (t,p ) = isinp n h(t,p), if 
P 7^ 0, the interface height become^ 

h(t,p) = h{0,p)e~ 2 ^-A= +J‘dr rjir,p)e- 2r (‘-^Mj^ (3.13) 
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The EW model is obtained by setting g(t) = 1 ++ %(t) = 0 [15] . 
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and the initial conditions (13. 5 p now read 

= H 0 S(p) , (h{Q,p)h{Q,q)^ = H 1 S{p + q) (3.14) 

The calculation of the long-time behaviour of physical observables will be based on eqs. (13. 1 1[) . (13.13jl . 
respectively. To do this in practice, g(t) must be found from the constraint (13.101) . with the 
explicit form (13. 1 21) of fit) taken into account. 


3.2 Linear responses 


Since the evolution of the model starts from a non-stationary initial state and the dynamics 
contains an external force which does not derive from a Hamiltonian, it leads to non-equilibrium 
relaxation, where responses and correlators must be analysed separately. We begin with the 
response functions. 

As the Arcetri model describes the growth of an interface, the conjugate variable to the 
interface height h is an additional particle-deposition rate j = J[s, q) at time s and at momentum 
q. Such an extra term can simply be added to the equations of motions by formally replacing 
rj i—>• rj + j. We then define the linear response of the interface height 


R(t,s;p,q ) : = 


6{h(t,P)) 


Sj[s, -q) 


= 0(f — s)5(p — q)‘ 


3=0 


d( s ) „- 2 r (t- 

9(t) 


2T(t—s)uj(p) 


(3.15) 


and used (13.131) . For a linear response as considered here, g(t) is still given by (13.101) . The 
Heaviside function 0(f — s) expresses the causality condition t > s. In direct space, this gives 


R(t, s;r — r') = 


5{h(t,r )) 


5j(s,r') 


3=0 


Q(t-s) I g{s) 
(2vr) d V gW 


dp e ip - 


ip-(r-r')-2T(t-s)u(p) = 0(t - S ) W F r-~ s) (3. 

V 9{t) 


16) 


where here and in what follows we use the abbreviation (with the modified Bessel function I n 
taken from [SSI eq. (9.6.19)]) 


Fr M : = 


(2rr) 


dp cos (p ■ r) e 


— 2 Truj{p) _ 


-n- 

a= 1 


- 2 rT 




(3.17) 


In order to derive the response of the slopes, consider first a source term in an action, 
expressed in momentum space 


d dp h(t, p)j{t, -p) 


d pJ2h(t, P )^j(t,-p) 

^ 1 Sill '[)„ 


a =1 


Sill p a 


d n 

/ d pu a (t,p)J a (t,-p) 

a= 1 d 


using (13.41) and where we introduced the integrated source J a (t,p ) := (— isinp a ) 1 ]{t,p). Con¬ 
sider the following linear response 


d 

Q{t,s\p,q) ■= Y 

a= 1 


S{Ug(t,p)) 
5Ja{s, -q) J=0 


6{p-q)Q(t- 



Y^ 2 Pa e-^d-sMp) 
a= 1 


(3.18) 
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which in direct space leads to 

Q(t,s-,r- r') = —^ [ dpdq e ip (r_r,) Q(t, s,p, q) 

( 2 tt y d J B2 

= / 8 d P A (p> elp ■ ( '' ^ ' , ■ 2^<, ■' ) “ (p, < 3 - 19 > 

In particular, the autoresponse of the slopes takes the simple form Q(t,s ) = Q(t, s; 0) = 
©0 - s)f((t - s)/2)yjg(s)/g(t). 


3.3 Correlation functions 


In order to write down the fluctuations and correlations of the height, we average the solution 
(I3.13P for p ^ 0. Then 

%p)) = (h(0,p)} acp( ^ (p)) (3.20) 

Going back to discrete momenta formulation for clarity, we have the deviation in the height 


h(t,r) - (h(t,r)) = 


N d 


E 


p i 


h{0,p ) - (h(0,p)) ) exp ( ^p ■ r - 2 Ttw(p)^ -^== 


27ri 


dr fj(r , p) exp ( — p ■ r - 2 T(t - r)w(p) 


' g(^) 

s(*) 


(3.21) 


such that the two-time height-height correlator becomes 

C(t, s;r — r') = C(s, t; r' - r) = (( h(t , r) - (/i(t, r)}) ( h(s , r) - (h(s, r')))) 

h( 0 ,p)/r( 0 ,-p)) - (h( 0 ,p) 


i_!_y 


^2mN 1 p-(r-r') g -2r (t+s)u(p) 


2 rr 

’"/v 5 " 


E 


r»min(t,s) 


dr 




\/9(t)9(s) 


^27rijV 1 p-(r*—r') 2r(£+s— 2t)uj(p) 


(3.22) 


For a vanishing initial width, see eq. (13.5ft with H\ = Fq , the term in the square brackets in the 
2 nd line vanishes for p = 0. We can now take again the N —> oo limit. Using F r (t) as defined 
in eq. (13.171) . the two-time correlator reads (simplified for t > s) 

h x 2rr 


C(t,s\r) = 


F r (t + s) + 


dr g{r)F r (t + s — 2 r) 


(3.23) 


VMM VWW) Jo 

and the autocorrelator becomes C(t, s ) = C(t, s: 0). In particular, the interface width reads 

HiF 0 (2t) 


w 2 (t) = C(t, t ) = 


2rr [ dr^F 0 (2t-2r) 
Jo 


9(t) Jo 9(t) 

Similarly, the slope-slope correlator can be found, again with t > s assumed 


(u a (t,r)u a (s,r')) = 




\/g{t)g{s ) W d J B 


dp sin 2 p a e 
9( T ) 

( 27r ) d Jb ^ Jo V g(t)g(s) 


ip-(r—r')—2T(t-\-s)uj(p) 


(3.24) 


(3.25) 


2rr 


dp / dr 


sin 2 p a e ip-b-r')-2r(t+ s -27>(p) 
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whereas ( u a Ub ) = 0 if a ^ b. We shall be interested in particular in the autocorrelator of the 
slopes, with f(t) given by eq. (I3.12p 


d 

s ) := {u a (t, r)u a (s, r)) 

a= 1 


^ /((^ + g )/ 2 ) 
1 y/9(t)g(s) 


+ 2TT [ dr 
Jo 


9(t) 

y/9(t)9(s) 


f 


Here, and in what follows, we shall always choose units such that T = 1. 



(3.26) 


3.4 Equivalence with the spherical spin glass 

As we now show, the ID Arcetri model is equivalent to the kinetic p = 2 spherical spin glass, 
analysed by Cugliandolo and Dean [64] 1^1 

The spin hamiltonian is H = Ji,j s i s ji where s* G M are spherical spins which 

obey the spherical constraint YliLi s f = N. The elements of the symmetric matrix J are 
independent gaussian random variables with zero mean and variance proportional to 1/N. Then 
the thermodynamic limit is well-defined and the probability distribution of the eigenvalues /i 
of J is given by the Wigner semi-circle law 

PM = ^ — (3.27) 

The corresponding eigenvectors of a spin configuration s(t) = {s*(t)} are denoted by s Al (t) = 
/i • s(t). Consider uniform initial conditions s^(0) = 1. The time-evolution is given by the 
Langevin equation [61, eq. (2.2)] 

dts^it) = (p + ^(t)) s„(t) + h„(t) + ^(t) (3.28) 

where h M (t) is an external magnetic held and ^(t) is the centered thermal noise, with variance 
(£n(t)£,u(t')) = 2TsG^5(t — t'). The solution of the Langevin equation is now immediate. 
Defining 7 (t) := exp(—2 f^dr^(r)), the spherical constraint can be rewritten as a Volterra 
integral equation [ 6 l, eq. (2.7)] 

70) = ((s/i( 0 ) exp (2 pt))) + 2 T sg [ dr q(r) ((exp 2 p(t - r))) 

Jo 

where the following average over the Wigner distribution p(p) is carried out [M] 

((exp( 2 /it))) := J dp p(n)e 2flt = ^ J dp. \/l - p 2 e 4/it = 

where [89] eq. (9.6.18)] was used. 

In order to see the relationship of these results with the ID Arcetri model, we first observe 
that the eigenvalue spectrum in the spherical spin glass is in the interval p G [—2,2], whereas 
the dispersion relation of the Arcetri model 2c a(p) G [0,4]. This can be matched through the 

14 A careful analysis [731190] shows that the ageing properties of truly glassy systems on one hand, and simple 
spin systems and the p = 2 spherical spin glass on the other hand, are quite distinct. 


(3.29) 


(3.30) 
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mapping /j 4 —2 + /i and 3 (f) 1 —>■ 3 (t) + 2 in order to keep the equation of motion (3.28) 
unchanged. Identifying the relationship 


g (t) = e 44 y(f) (3.31) 

it follows that the Volterra equation \3.29i) with the explicit average \3.3C\) becomes exactly the 
spherical constraint H3.10 \ ) . where f{t) is given in (13.12(1 and the identifications 


T = 2T sg and Hi = 2 


(3.32) 


Next, the magnetic autoresponse of the spin glass (with t > s assumed) ( 6 H eqs. (2.16,3.17)] 


N 

Rsg( t, s ) := 'Y 

i —1 


5hi(s) 


h =0 



((exp n(t - s))) 



2 Q(t,s) 


(3.33) 

is identical to the slope autoresponse (I3.19jl . up to a factor 2. Finally, and recalling (13. 3 Oh the 
disorder-averaged spin glass auto correlator [64, eqs. (2.12,3.6)] 


Csg(M) 


N 


Y ( Si(t)Si(s )) 


2—1 


dpp(/i) (s M (t)s M (s)) 


(3.34) 


' —2 


ViWKi 


((exp pit + s))) + 2 T sg 


^min(t,s) 


dr y(r) ((exp ji(t + s - 2 r))) = 4(t, s) 


reduces to the slope-slope autocorrelator (I3.26p . with the same identihcations (I3.32p as above. 

Since it is well-known that the spherical spin glass is in the same universality class as 
the 3 D kinetic spherical model [64j [65], it follows that the ID Arcetri model is in the same 
universality class as well. The physical correspondence is between the slopes in the Arcetri 
model and the magnetic spherical spins, see (I3.33p . (l3.34p . This can be generalised: the long¬ 
time asymptotics of g{t ) and fit) imply that the exponents which give the long-time behaviour 
of Aft, s ) and Q{t, s) of the d-dimensional Arcetri model are the same as those of the magnetic 
autocorrelators and autoresponses in the kinetic spherical model 165\/ in d + 2 dimensions. 


4 Long-time behaviour 


Given the form of the spherical constraint (13.101) as a Volterra integral equation, its solution 
proceeds via standard Laplace transforms [SXUMi .53- Let ~g{p) = C{g{t))[p) = dt g(t)e~ pt 
denote the Laplace transform. Then (13. lOj) can be inverted 


ff(p) 


Hifjp) 

d~2Tf(p) 


(4.1) 


Standard Tauberian theorems [1091 ch. XIII.5] relate the long-time behaviour of g(t) to the 
p —> 0 behaviour of ~g{p). Hence, expanding f(p) for p 1, one can find first g(p) and then 
g{t ), for t large enough. The explicit calculations are detailed in appendix A. 
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4.1 Fast relaxation for T > T c 


Consider the situation when the denominator in (14.1 ft vanishes for some value po > 0: d — 
2Tf(po) = 0. Since f(p) decreases monotonously with p, this zero exists and is simple for T 
large enough. Hence ~g{p) has a simple pole at p = pq. Transforming back, one has 


g(t) = 


Hid 


aP0t 


2 T f'(p o)' 4 T 2 f'(p o) 


(4.2) 


Here l/po is the finite relaxation time and governs the time-translation-invariant approach 
towards the stationary state (tco,i are constants) 


~ r d/2 exp 


' 2 P0T " 47 


(4.3) 


R(s + r, s; r) 
w\t) 

C(s + r, s) 

This rapid relaxation does not show dynamical scaling and is of no particular interest to us. 


w 0 T + w\e pot 

1 


exp 1 -- Pot 


4.2 Ageing at the critical point 


The smallest temperature T = T c , for which g(p) has a singularity, occurs when d/ (2 T c ) = /(0). 
Hence the critical temperature T c = T c (d) is given by 

1 2 _ 1 r°° 

m = -^ )=2 L dte (4 - 4) 

Clearly, T c (d) > 0 for all d > 0. Using [89, eq. (11.4.13)] and [HI eq. (2.15.20.6)], we have the 
explicit valued 

2n 

T c ( 1) = 2 , T c { 2) =-- ~ 5.5038 ... , T c (3) ~ 9.53099 ... (4.5) 

71 — 2 

In addition, we have T c (d) ~ d for d 1 and T c (d) ~ 4 d for d 1. 

We now write down the linear response (I3.16|) . the interface width (I3.24j) and the autocor¬ 
relator (I3.23p . in the long-time scaling limit, where both t, s —> oo such that y = t/s is kept 
fixed. We begin with the two-time time-space response and expand g(t) and F r (t) for large 
times. This gives 


R(t,s\r ) = 

R(t, s ) = 
with the explicit scaling function 

fn{y) = (4t r )- rf/2 (z/ - l)- rf/2 2 /( 2 -^)/ 4 ; if 0 < d < 2 

fR(y) = (47r)” d/2 (y - l)~ d/2 ; if 2 < d 

15 For d = 1 in agreement with the spin glass critical point Xsg, c = 1 5T - 


'g{s) 


g(t) 

s~ d/2 f R (t/s) 


F r (t — s) ~ R(t, s ) exp 


1 r 2 


At — s 


(4.6) 


(4.7) 
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From this, the values of the non-equilibrium exponents a, \r and also the dynamical exponent 
z — 2 can be read off and are listed in table [H 

Next, consider the two-time correlator. In the scaling limit, with t = ys > s, and for 
0 < d, < 2, we find 


C(t, s; r) ~ 2Tdts)^ 2 d / dr/’ 1 F r (t + s — 2r) 

Jo 

s 1 -"/ 2 27) y^ 4 t. 4 ' 2 - 1 (4tt [(y + 1) - 2u])-“ /2 exp 


(4.8) 


s y + 1 — 2u 


which is of the expected scaling form C(f,s;r) = s' -’’FcO/sy/s). 
calculation gives 

C(t,s-,r) = s 1 ~ d/2 2T c [ du u d/2 ~ l (47 t [{y + 1) — 2u])~ d/2 exp - 

Jo 


For d > 2, an analogous 


s y + 1 — 2u 


(4.9) 


Especially, the autocorrelator C(t,s) = C(f,s; 0) = s b Fc(t, s; 0) = s b fc(t/s ) has the explicit 
scaling function^ 


fc(y) 

fc(y) 


4 T c {d) t/ (2 ~ d)/4 

d(47r) d / 2 (1 + y) d / 2 

4 T c (d) 

d(47r) d / 2 




2 2 


2’1 + y 


y 


— (3d—2)/4 


(! + J/) 


—d/2 


F , d d ., , d . 2 
2jFi1 2'2' 1+ 2'TT^ 


—d/2 


if 0 < d < 2 
if 2 < d (4.10) 


along with the asymptotics for y large. The values of the exponents b and Ac are listed in 
table [0 

In analogy with the known results of the kinetic spherical model of magnets HUGS], the 
behaviour at the upper critical dimension d = d* = 2 of the two-time observables R(t, s; r) and 
C(t,s\r) can be obtained by formally taking the limit d 2 in eqs. (14.71) and (14.101) . Hence 
no logarithmic modifications of the leading scaling behaviour arise, but there are additive 
logarithmic corrections to scaling. This is different for single-time quantities like the interface 
width w(t), of which the leading long-time behaviour is 


w\t) 

27 xT c {d) O —d/2 . 

if 0 < d < 2 

(87r) rf / 2 sin(7rd/2) 


w\t) 

~ 2 T c In t ; 

i-b 

Si- 

to 

w\t) 

= 2T c (d)F d + 0(t 1 ~ d ' 2 ) ; 

if d > 2 


with the constant F ( j := / 0 °°dr (e 4r /o(4r)) d , which diverges for d \ 2@ Hence, the inter¬ 
face width w(t) ~ Vint growths logarithmically in 2D , see table CD in agreement also with 
experimental findings [40). The values of exponent (3, listed in table [T] are the same as for the 
EW-universality class. 

Inspection of the values of the various exponents collected in table Q] shows that 

16 Recast the integral fjdu u d / 2_1 \{y + 1) — 2 u]~ d ^ 2 = |(1 + y)~ d F 2 F\ ( d/2, d/2 ; 1 + d/2 ; 2/(1 + y)) 

= (y — l)~ d t 2 B 2 /(\+ y ) ( d/2 ,1 — d/2 ) as a hypergeometric or incomplete Beta function [551 eqs.(15.3.1,26.5.23)]. 

17 T 3 = n~ 2 (18 + 12-v^"— lOVT— 7^/6") K 2 [(2 — v / 3)(>/3 — >/2)] ^ 0.1263655..., where K(fc) is the com¬ 
plete elliptic integral [91] eq.(2.15.21.2)]. 
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1 . the stationary exponents z, (3, and consequently also the roughness exponent a = z/3 , 
of the Arcetri model are in all dimensions the same as those of the Edwards-Wilkinson 
model. This might have been anticipated, since the equation of motion (13.2[) is still linear 
- and also reflects properties of the magnetic spherical model, where the equilibrium 
critical exponent /3 — \ of the order parameter and 77 = 0 of the spin-spin correlator keep 
their mean-field values [59] . 

2 . the non-equilibrium autocorrelation/autoresponse exponents Ac = Xr depend on the 
dimensionality in a non-trivial way. Only for d > d* = 2, their values are the same as in 
the EW-universality class. 

In this respect, the dynamics of the Arcetri model is quite similar to what is found in the 
non-equilibrium critical dynamics of simple magnets with a non-conserved order parame¬ 
ter and disordered initial states. For such systems, a long-standing result of field-theory 
[13 , 21} IT] asserts that the true non-equilibrium exponents Ac, Xr should be independent 
of the stationary exponents (since an independent renormalisation is required for their 
calculation |13|). The role of the stationary exponents is taken here by z,a,(3 and the 
exponents a, b related to them. Hence the Arcetri model with d < 2 can be considered as 
an explicit example of this general fact. 

We also recall that in the kpz class with d < 2, a Ward identity prevents this addi¬ 
tional renormalisation, to all orders in perturbation theory [ 22 ], which is a qualitatively 
different situation. We are not aware of any estimate of Ac from a non-perturbative 
renormalisation-group study in the kpz universality class, for d > 2 . 

3. for the KPZ-universality class, there is a conjecture^ for a flat interface, that A c = d 
[IT] ; this also happens to be satisfied in the EW-model. However, this conjecture does not 
extend to the Arcetri model, since Ac = Xr d for d < 2. 

4. in the Arcetri model, autocorrelation and autoresponse exponents are always equal. Is 
the available numerical evidence for Ac f Xr in the 2 D KPZ-model [27] the final word ? 


Remarkably, even for d > 2, the Arcetri and the EW models are in different universality classes, 
in spite of all their exponents being equal. In order to see this, recall that because of a = b, one 
can define the fluctuation-dissipation ratio (FDR) [63] (for reviews, see [92, 2T] |93} -9!} 195] ) 

X(t,s):=TR(t,s)^AffAj =*(}) (4.12) 


where the last relation holds in the scaling limit. In magnetic systems or spin glasses, one uses 
the value of X(y ) — 1 as a measure of the distance with the respect to an equilibrium state. 
According to the Godreche-Luck conjecture [65], the limit FDR X^ = lirn^oo X(y) should be 
an universal number. For the critical Arcetri model, the explicit results (14.711 . (14.10p lead to 


f d/(d + 2) ; if 0 < d < 2 

\ d/4 ; if d > 2 


(4.13) 


18 For the kpz class in d < 2 dimensions, Ac 


d has been proven to all orders in perturbation theory [22] . 
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Figure 2: Fluctuation-dissipation ratio X(y) over against y — t/s of the critical Arcetri model 
for (a) d < 2 and (b) d > 2. The small circles on the right axis indicate the limit ratio X^, eq. 
114. 1311 . In (b), the solid green line gives X(y) for the EW-model in d — 4 dimensions. 

Eq. (14.13ft is distinct from the well-known value ^ of the EW-model, valid in any 

dimension [15] 0 To illustrate this further, we show in figure [2] the scaling function X{y) as a 
function of y = t/s, for d below (figure^) and above (figure [2]o) the upper critical dimension 
d* = 2. Starting from the value AT(1) = 1 at equal times, the limit ratio X is continuously 
approached for an increasing temporal separation y = t/s —> oo. Qualitatively, this approach 
is monotonous for dimensions d < 2, but for d > 2 but not too large, X(y) goes through a 
maximum before approaching its limit value. For sufficiently small dimensions, this approach 
is from above, but when d becomes large enough, the maximum of X(y) disappears and the 
limit value X^ is approached monotonously from below. Curiously, for d — 4 one has X <*> = 1 . 

For d = 2, the FDR of the EW-model is identical to the Arcetri model and the green solid 
line in figure [ 2)3 gives the FDR of the EW-model for d — 4. For all dimensions 2 < d < 4, the 
FDRs of the EW-model fall between these two curves, which are quite distinct from those of 
the Arcetri model. 

The distinction between the Arcetri and EW-models is further illustrated in figure [3) where 
the ratio of the two fluctuation-dissipation ratios is plotted. With an increasing separation 
y = t/s > lof the time-scales the distinction of the two models becomes progressively more 
clear. Only for d = 2 the two functions X(y) are identical, hence the two models are indeed 
identical in 2 

Hence the Arcetri model with d > 2 illustrates the usefulness of the fluctuation-dissipation 

19 The slope autocorrelator A(t, s) and autoresponse Q(t, s) of the d-dimensional Arcetri model have the same 
scaling behaviour as the magnetic autocorrelator and autoresponse in the (d + 2)-dimensional spherical model, 
see sec. 3.4. For d < 2, the limit ‘slope FDR’ A™(d + 2) = 1 — 2/(d+ 2) = d/(d + 2) [65] agrees with the height 
FDR (14.1311 of the d-dimensional Arcetri model. 

20 Therefore, the experiments described in [40] apply to both. 
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Figure 3: Ratio of the fluctuation-dissipation ratio X\ TC (y) for the Arcetri model and of Xew(u) 
for the EW-model, as a function of y = t/s, for several dimensions 2 < d < 4. 

ratio as a diagnostic tool for a fine distinction of non-equilibrium universality classes. It also 
points to the interest of measuring experimentally both two-time correlators and responses in 
growing interfaces (long-range interactions generically reduce the value of d*). 

4.3 Relationship with the bosonic pair-contact process 

Surprisingly, it turns out that for 2 < d < 4, the non-equilibrium relaxation properties of 
the critical Arcetri model are related to those of a different stochastic process, the so-called 
(multi-)critical bosonic pair-contact process with diffusion (BPCPD) [IFIl 1971155] . see p3J for a 
field-theoretic description. In this model, each site of a d -dimensional hypercubic lattice can be 
occupied by an arbitrary number n G N of particles of a single species A. Single particles may 
hop to a nearest-neighbour site, with a diffusion rate D. On the same site, pairs of particles 
may react, according to the schemes 2A —>■ (2 + k)A or 2A —> (2 — k)A, where k is either one or 
two. The universal long-time behaviour of the model does not depend on the value of k. The 
two reaction rates are chosen equal T[2A —>■ (2 + k)A] = Y(2A —> (2 — k)A] = /i. The model’s 
behaviour is determined by the value of the control parameter a := k 2 /i/D. There is a critical 
value, a = ac, witb 0 

1 f 00 

— = 2 d u (e~ Au UAu)) d (4.14) 

ac Jo 

and for d > 2, one has ac > 0. For a > ac, all particles cluster on a few sites only [97J [98l E| 
while for a < ac, the system evolves towards a spatially homogeneous stationary state. We are 
interested in the multi-critical point at a = ac > 0, where this ‘ clustering transition’’ occurs. 
For d < 2, the multi-critical point in the BPCPD does not exist. 

21 Using from (14.111) . ac = ( 2 .F 3) -1 ~ 3.956776... for d = 3. 
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The dynamics of the model is described in terms of the bosonic particle-annihilation operator 
a(t, r ) such that the particle-number operator is n(t, r) = a'(t, r)a(t , r). For any value of a, and 
an initial state with uncorrelated particles of a mean number density po, the average particle- 
number (■ n(t,r )) = ( a(t,r )) = po is constant [971 [98] . The two-time density-density correlator 
can be expressed as [9] 


(n(t, 0)n(s, r)) - pjj = C(t, s; r) + R(t, s; r)p 0 
with the following definition of the two-time correlators and responses 

8{a(t, 0 )) 


C(t,s;r) = (a f (f,0)a(s,r)) - pi , R(t,s;r) = 


Sj(s,r) 


(4.15) 


(4.16) 


i=o 


where j denotes the particle-creation rate conjugate to the particle number. The usual scaling 
behaviour (ll.2p . fll.3p has been confirmed, for all a < ac- Precisely at the clustering transition 
a = and for 2 < d < 4 dimensions, the explicit scaling behaviour reads, for s —> oo [98] 


C(ys,s;0) = c 0 pl(y+ 1) d/2 2 F l (4 ll + 4 |2_A , R(ya, s; 0) = s d/2 r 0 (y - 1) !i/2 

(4.17) 

and where Cq, r 0 are known normalisation constants. Hence, C(t, s; 0) can be considered as the 
scaling limit of the connected density-density autocorrelator in the BPCPD. For 2 < d < 4, 
the shape of the scaling functions is identical to the ones of the Arcetri model given in eqs. 
(I4.10p . (l4.7p . However, although a — d/2 — 1 in both models, the values of the ageing exponents 
are distinct, namely 6 = 0 for the bpcpd and b = d/2 — 1 for the Arcetri model. Since a ^ 6, 
the fluctuation-dissipation ratio (14.121) cannot be defined in the multi-critical bpcpd. 

This is a further example which illustrates that the values of the non-equilibrium exponents 
A c,r are independent from the stationary exponents a,/3,z, analogously to non-equilibrium 
field theory of critical magnets pGSEUIZ!, but different from the kpz universality class |22l 13- 
However, this independence is illustrated in the Arcetri model for 2 < d < 4 in a new way. 
Recall that for d < 2, all the exponents a, (3, z and a, b are found to be the same in the EW and 
Arcetri universality classes, yet the values of A c = A r are different. Here, when 2 < d < 4, 
the values of the non-equilibrium exponents Ac = \r = d, as as well the ageing exponent 
a — d/ 2 — 1, are the same. However, the ageing exponent b = 0 in the bpcpd and b = d/ 2 — 1 
in the Arcetri model are different. 


Finally, for d > 4, the shape of fc{y) hr the multi-critical BPCPD is different from (14.17p . 
although Ac = d, and b = d/ 2 — 2 p5]. This gives a different example, but of the same kind, 
as found for d < 4. 


We emphasise that although the value of the autocorrelation exponent Ac = d, the two-time 
correlation functions are distinct from those of the EW-model, where C(t,s ) = s~ b f^ W \t/s) 
with b = d/2 — 1 and f^ W \y) = ^^(47r) _d//2 [(y + l) 1-d / 2 — | y — l] 1 -^ 2 ] [15]. This rather 
coincides with the shape of fc{y) in the critical bpcpd with a < ac [98]. 


4.4 Ageing for T <T C 

In order to discuss the behaviour below the critical point, when T < T c and the noise is thus 
relatively weak, it is useful to introduce the auxiliary variable m 2 := 1 — T/T c . 
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The two-time response function again takes the scaling form (14.61) . where the scaling function 
fn(y) of the autoresponse now becomes 

h(y) = (‘tx)- d, V‘ +2>,i (y - ir d ' 2 (4.18) 

for all dimensions d > 0. Similarly, the two-time autocorrelator becomes (with t = ys > s 
assumed; the first correction is valid for d < 2 , see appendix A) 

C(t, s ) = 2 d+2 m 2 sy {2+d)/A (y + l)~ d/2 + O (! Ts l ~ d/ 2 ) (4.19) 

which has the expected scaling form C(t, s ) = s~ b fc(t/s ) with the explicit scaling function 

fc(y) = 2 d+2 m 2 y {2+d)/4 (y + l)~ d/2 (4.20) 

The non-equilibrium exponents can be read off and are listed in table [TJ The extension to time- 
space correlations C(t,s;r) is straightforward. The interface width is obtained as the special 
case w 2 (t ) = C(t,t) and reads (the first correction is valid for d< 2 ) 

w 2 (t ) = 4 m 2 t + O (Tt ] - d/2 ) (4.21) 

For correlators and the interface width, the long-time behaviour only depends on the ratio 
T/T c and the short-ranged nature of the initial correlations (parametrised by Hi). The ‘thermal’ 
convolution term, which dominated the critical case, merely gives rise to a correction to scaling. 
This leading result w 2 {t ) ~ 4 m 2 t, valid for t large enough, is independent of d, hence the growth 
exponent /3 — |. 


4.5 Global persistence 


A different class of observables are first-passage probabilities, such as persistence probabilities, 
see [ 100] for a recent detailed review. The global persistence probability Pcit) is the probability 
that a fluctuation 5h G (t ) of the global height h G (t) f Rd dr h(t, r ) h(t, 0) did never change 
its sign between the initial instant t = 0 and time t. For large times, one expects a power-law 
decay Pq( t) ~ t~ dG , where 6 G is the global persistence exponent. In the infinite-size limit, the 
central limit theorem implies that the global height should be a gaussian variable I 111 1 . 100 ]. 
Remarkably, it can be shown that then the model’s underlying stochastic process is markovian 
if and only if the normalised global autocorrelator 


N(t,s) 


C 0 (t,s ) 


C 0 (t,t)C 0 {s,s ) 


- if)' 


(4.22) 


(where Co(t, s) is the two-time autocorrelator in momentum space) takes in the scaling limit, 
t, s —* oo with y = t/s fixed, a simple power-law form jlOlj . Then 0 G = /i, and in addition, a 
scaling relation relating 0 G and the autocorrelation exponent Ac can be derived, for quenches 
either to T = T c [101] or to T < T c [1021 103 . ilOOj . Adapting this to the exponent notation 
(II.21) . (II.31) used for growing interfaces, this ‘markovian’ scaling relation reads 


0 G 




(4.23) 
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where the bound follows from a Yeung-Rao-Desai inequality [ I'iTdj . see appendix B. However, 
in many non-equilibrium systems, quenched to either T = T c or T < T c , those scaling relations 
turn out to be broken (albeit by numerically small amounts) so that one should conclude that 
the underlying stochastic processes cannot be markovian at sufficiently large times, see !H1 
and references therein for explicit examples. Practically, deviations of N(t, s ) from a pure power 
law may be more easy and more reliable to detect than tiny deviations from (14.231) [103] . 

Applying this general method to the Arcetri model, with (13.131) the global height autocorre¬ 
lator is Co(t, s ) = [Hi + 2rr J 0 s dr g(r)) / y/g(t)g(s) , for t > s. Recalling the leading long-time 
behaviour of g(t ) from appendix A, we have 


H i + 2VT Ip drg(r) oo /s \ ** 

H 1 + 2rTf 0 t drg(r ) KtJ ’ M G 

As expected, the verification of the simple power-law 
hence the identification g = 9 q is admissible. 


N{t,s) = 


{ 0 

; if T <T C and d > 0 

= { d/4 

; if T = T c and 0 < d < 2 

1 1/2 

; if T = T c and d > 2 


(4.24) 

(|4.22|) confirms the Markov propertv. 


5 Discussion and perspectives 

Concluendum est: this work explores the properties of an analogue of the familiar spherical 
model of ferromagnets [59] for the description of growing interfaces. In the scheme presented 
here, we admitted that the correspondence should be made between the local slopes of the 
interface and the spherical spins, ft then becomes relatively straightforward to write down and 
solve the Langevin equations of motion, since essentially all techniques can be borrowed from 
the kinetic spherical model pH, 64||65]. In particular, we have analysed the long-time behaviour 
of the Arcetri model, defined in section 2. Our findings are as follows: 

1. the parameter T, named a ‘temperature’ by analogy with the kinetic spherical model, 
and defined through the noise correlator (12.21) . admits for all dimensions d > 0 a critical 
value T c = T c (d), given by eq. (14.4|) . Qualitatively, the Arcetri model is in this respect 
more analogous to magnetic systems and qualitatively different from more usual interface 
growth universality classes, such as described by the kpz or EW universality classes. 

At criticality T = T c [d), starting from a flat substrate, the interface becomes rough for 
d < 2 and remains smooth for d > 2. For T < T c [d), the interface always becomes rough 
and for T > T c [d), even an initially rough interface becomes smooth in a finite time. 

2. for T > T c [d ), the model’s relaxation behaviour is rapid, time-translation-invariant and 
governed by a finite relaxation time. 

3. in contrast, for T < T c , the relaxation is slow, time-translation-invariance is broken and 
the two-time observables obey dynamical scaling (jl.2j) . (jl.3[l for large times. Hence the 
conditions for physical ageing, as defined in [9], are obeyed. 

Such a physical ageing behaviour has already been found in more standard models of 
growing interfaces, especially in the KPZ- and EW-models. In table [D the values of the 
relevant exponents are listed. 
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4. the critical Arcetri model has an upper critical dimension d* = 2. 

If d — d* = 2, no logarithmic modifications of the dynamical scaling arise for the two- 
time observables, but a logarithmic behaviour is seen for single-time quantities such as 
the interface width. The 2D Arcetri and EW-models are identical. 

5. in many respects, if one chooses T = T c , the resulting behaviour of the critical Arcetri 
model is qualitatively quite analogous to what has been found in other universality classes. 
However, some differences can be seen as well: 

(a) the stationary exponents z,/3,a (as well as the exponents a,b which depend on 
them), are the same for the EW-model and the Arcetri model. However, for d < 2, 
the non-equilibrium exponents Ac and Xr of both models are different. 

In this respect, the critical Arcetri model behaves analogously to the non-equilibrium 
critical dynamics of simple non-conserved magnets and is qualitatively different from 
growth models in the kpz universality class [22]. It provides an explicit example 
where, for disordered initial states, the non-equilibrium exponents A c,r are indepen¬ 
dent (i.e. not related by a scaling relation) from the stationary exponents. This 
is a long-standing prediction from the non-equilibrium critical dynamics of simple 
non-conserved magnets mmm- 

(b) in the Arcetri model, one has always Ac = Xr. This matches the exact results of the 
EW-model and the ID KPZ-model. Curiously, the available numerical data indicate 
the contrary for the 2D KPZ-model [27]. 

(c) the conjecture A c = d \22, E], formulated for the kpz equation, and observed to 
hold in the EW-model as well, does not extend to the Arcetri model when d <2. 

The available numerical data for the 2D KPZ-model [271 EE] do not agree with it 
either. We remark that for the ID KPZ equation, the shape of the two-time re¬ 
sponse function could only be explained using the logarithmic extension of local 
scale-invariance, which implies logarithmic terms in the scaling function [181 [105] . If 
such an observation could be extended to the 2D case as well: could un-recognised 
logarithmic contributions have modified the effective values of Ac or A r ? 

6. for dimensions 2 < d < 4, the two-time response function, and the shape of the scaling 
functions of the two-time correlator (up to normalisation), are the same as at the multi- 
critical point in the so-called bosonic pair-contact process with diffusion (bpcpd) [98j . 
However, although the asymptotic exponent has the same value Ac = d, the shape of the 
scaling function of the two-time correlator of the EW-model is distinct from the one of 
the Arcetri model for d >2, although the two-time responses R(t, s; r ) and the exponent 
b = d/2 — 1 agree. 

7. an universal method to distinguish the Arcetri and EW-models for d > 2 uses the 
fluctuation-dissipation ratio (14 .12 j) [63] . in particular the universal limit FDR X [65] 
has different values in the two universality classes J^l 

Hence, for d > 2, we have three distinct models, with the same values of the exponents 
a, and A r = Ac, but which are distinguished as follows: (i) between the Arcetri and 

22 The distinction between the two models for d, f 2 should be related to the spherical constraint (12.41) . which 
is incompatible with a flat surface h n = cste.. 
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bpcpd- classes, the exponent b has different values, (ii) between the Arcetri and EW- 
classes, the FDR and in particular the limit are different. This illustrates once more 
the independence of the exponent Ac from the other ones. 

8 . the relationship with the kinetic spherical model in d + 2 dimensions (see section 3.4), 
relating local slopes to local magnetisations, allows to draw a clear physical picture on how 
the long-time relaxation is going on: for T < T c , the Arcetri model undergoes a coarsening, 
such that increasingly larger patches of the interface (of typical size L(t) f 1 / 2 for all d 
and T < T c [61) [79]) have constant slope, leading to a saw-tooth pattern. For T = T c , 
the interface slopes are correlated over increasing distances L(f), such that the interface 
itself should become a fractal. 

The spherical model is easily adapted and generalised, for instance to long-range initial 
conditions [69] EE], long-range interactions (S3 E51 EB], external fields im conserved order 
parameters [621(631711177!, frustrations and/or disorder [64], external drives (BU [H2] and so on. 
Several of those modifications could be of interest to experimentalists. It should be possible 
to generalise the Arcetri model in these directions, to work out the consequences on its ageing 
behaviour and explore relationships with long-range particle-reaction models [1116] . 

Another important question will be if the Arcetri model can be understood as a n —> oo 
limit of a suitable n-component generalisation of the KPZ-equation, following the lines outlined 
hi [ 1 TT 7 ] . We hope to come back to this elsewhere. 

Finally, the model analysed here is but one way to write down a ‘spherical analogue’ of the 
kpz universality class. We used here the representation in terms of the Burgers equation 

d t u = vd 2 x u + /md x u + d x r] H > d t u = vd 2 x u + i(t)u + d x r] (5.1) 

where 3 (f) (dx'u) 1/t (which follows from g(t) ~ t F ) might be seen as some kind of 
‘averaged curvature’ of the interface. However, the choice considered here was also motivated 
by its mathematical simplicity. It might be thought that a more faithful representation of the 
structure of the non-linear terms could be achieved in terms of a replacement 

d t u = vd x u + fiud x u + d x r) *->■ d t u = vd 2 u + i{t)d x ii + d x 77 (5.2) 

where 3 (f) ~ (u) might be viewed as some kind of ‘averaged slope’. I 11 these two cases, the 
Lagrange multiplier 3 (t) is to be found from a spherical constraint ^2 x {u 2 ) = J\f. Finally, we 
might have started directly from the KPZ equation 

d t h = vd x h + (d x h) 2 T rj h-)- d t h — vd 2 h + l(t)d x h + rj (5.3) 

where 3 (t) ~ (d x h) might again be interpreted as an ‘averaged slope’ and will be found from 
a constraint ^2 x {(d x h ) 2 ) = A f. While eq. (15.11) has been studied here, work on the other two 
models, as defined in eqs. fl5.2H . fl5.3H is in progress. The analysis of their sphericl constraints 
presents new features not seen in the system studied in this work and apparently leads to 
new types of behaviour. A sequel paper will take up the analysis of these distinct universality 
classes. 
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Appendix A. Computation of the long-time behaviour 


The explicit long-time behaviour of the spherical constraint parameter g(t), for T <T C and the 
consequences for the interface width w(t) and the two-time responses and correlators, stated 
in section 4, will be derived. 

Since g(p) is given by (14.11) , we use standard Tauberian theorems [1091 ch. XIII.5], in 
order to obtain the long-time behaviour of g(t ) from the p —> 0 behaviour of ~g{p). To do so 
explicitly, two methods have been used: (i) one may consider g(t),f(t) as regular functions, 
which are found by inverse Laplace transformation. These must be completed by certain sum 
rules, needed for the computation of the correlators, see e.g. [63, ESI EZ1 EJj. (ii) perhaps 
more straightforwardly, accept that g(t),f(t) may have singular terms described by the Delta 
function and its derivatives, which are formally obtained by inverting the first few orders of 
g(p) [OUZO], see below. Although these singular terms do not contribute to the responses, they 
are important when inserted into the integrals which give the width or the correlators, since the 
singular terms in g{r) will reproduce the effect of the sum rules just mentioned. In conclusion, 
both methods lead to the same final result. 


A.l Ageing at the critical point 


Here, we consider the long-time behaviour at the critical point T = T c (d), as given by (14.41) , (14.51) . 

In order to find g(p), we first expand f(p) = J 0 °°dt e~ pt f(t), where f{t) was given in (13.12)1 . 
In principle, such an expansion will contain terms analytic in p, i.e. ~ p n with n G N, and 
non-analytic contributions, i.e. ~ p a with a N. The origin of both is clearly seen by splitting 
the integral J 0 °°d£ = f^dt + J^°dt where p is a cut-off to be sent to infinity at the end. The first 
term can be formally expanded in p, as long as the corresponding coefficients converge, and 
gives the analytic part. The non-analytic part is found by inserting the t —> oo asymptotics of 
fit) in the second term. To carry this out this, one sends first p —> 0 and only afterwards, one 
recovers the p —> oo limit. This standard calculation, e.g. [ 831 11081 169] . leads to the expansions, 
with A~ l = -f(8vr) rf / 2 sin (f) T (f), /(0) = d/( 2 T c (d)) and /'(0) = f™dte- dt I 1 (t)I 0 (t ) d - 1 


2 j^(d)+A d p d/2 ; if 0 < d < 2 

7(P)^{ 2 OT+nbOnP + C^- 1 ) i Hd = 2 
^ + /'(o)j) + V /2 


(A.l) 


if 2 < d < 4 


(Ce = 0.5772... is Euler’s constant) and for d > 4, a term O (p 2 ) appears which does not 
contribute to the leading scaling behaviour. Hence, with G d ■— — 7 > 0 

\ z± c{a)) f (0) 


9 ip) - 


(2 T c (d)) 2 A d P 


H 1 d/2 


Hi 

2 T c (d) 


+ 0 {p) 


f (0) 
if 0 < d < 2 


n , _ Mi _|_ H 1 A d _ 1_ d/2—1 1 . if 9 < cl < 4 

Ud P 2 T 0 {d) ' (2 T c (d))' 2 (/'(O)) 2 ^ + ’ ’ 11 z ^ a ^ ^ 


(A.2) 


and inverting this term-by-term gives, for t —> oo and with g d := — ^rld)) 2 A 7 r(!/ 2 ) > ® 


9(t) ^ 


9d 


t d / 2 


-1 


Hi 


m 


if 0 < d < 2 


G d 


2T a {d) 

+0{t~ d / 2 ) ; if d > 2 


(A.3) 


2 T c (d) 
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In order to find the two-time response explicitly, it is enough to insert into (j3.16jl the 
asymptotic form of g(t) for t —> oo, as derived above, to simply drop any distributional tenrrl 22 ! 
and to expand F r (t) for large times as well0 This immediately gives (14.61) . (14.711 . 

Next, we analyse the interface width w(t). Consider first the case 0 < d < 2. From (I3.24p . 
it follows 


w\t) = 


W) Fo(2t) + W) / dT » (T)jr « (2 ‘- 2T) 


W) Fo(2t) + l l/ dT 


Hi 

2T r 


- nsr */ 2 - 1 


S(t) - g d r' 


F 0 (2t - 2r) 


~ 2T c t 1 ~ d/2 / dr r d/2_1 F 0 (2t — 2r) 

Jo 

= 2T c t 1 ~ d/2 C~ l ((£r d / 2 - 1 ) (p) (CF 0 (2 t)) (p)) (t) 

2 t /( 1 ~ d/2 ^ d/2) ( C - 1 -) ( t ) tJ - d / 2 

Pj 


( 87 r ) d / 2 

2vr 


(87r) d / 2 sin(7rd/2) 


T c (d)t 


1—d/2 


(A.4) 


which gives the first line in (14.lip . In the second line in (lA.4p . we used in the integral the 
asymptotic form of g(r), including both the singular distributional as well as the leading regular 
term. Integrating the singular term, we see that it cancels against the other contribution to 
w 2 {t). In the third line, we also inserted the asymptotic form of g(t). As for the response before, 
the singular terms do not contribute for t —> oo. Furthermore, the non-universal amplitudes g d 
contained in g(r) and g{t) cancel. In the forth line, we recognise the remaining integral as a 
Laplace convolution which is treated via Laplace transforms. In the fifth line, we inserted the 
small-p behaviour of the two Laplace transforms, namely (.Ct^ 2 ^ 1 ) (p) = T(d/2)p ~ d / 2 and 


d P-s-0 


(CF 0 l2tj) ( p ) = / it e-* (<: '7 0 !4ti)“ 


d t e--‘(Snty^ + 0(1) ”3° 12L-A{2± pp-' 


In analogy with the computation above of f(p), the integral diverges for d < 2 in the p —> 0 limit. 
Hence its leading contribution in this limit can be found by splitting the integral J 0 °° = f d + f^° 
with a cut-off r). The contribution of the first term remains finite for p > 0 and rj < oo and 
it is enough here to extract the leading, divergent, behaviour of the second term. Herein, the 
T-functions are defined via analytic continuation, if necessary [89]. On the other hand, for 
d > 2, we first find the same cancellation of the ^-dependent term and then, using (1A.3P 

pt poo poo 

w 2 {t) ~ 2T C / dr F 0 (2t - 2 t) = 2T c / dr F 0 (2r) -2T C / dr F 0 (2r) = 2 T c T d + 0(t l ~ d/2 ) 

Jo Jo _ / Jt 

= -Td 

(A.5) 

since the amplitudes G d cancel. The last integral was estimated by using F 0 (t) ~ t d / 2 , for t 
large enough. This gives the last line in (14.lip . 

23 The singular terms in (IA.3II do not contribute for t large. Heuristically, recall: 8(t) = lim^^oo \/A/7J- e _A * 2 
vanishes for t ^ 0. 

24 The relevant Bessel asymptotic formula is I r (t) ~ ( 27 rf) _1 ' /2 e 4_r2 /l 2 b(l -g 0(t _1 ) for t, —> oo. 
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For the two-time correlator, consider first the case 0 < d < 2. We can use the same can¬ 
cellation mechanism as before for the temperature-independent term and then find, assuming 
t = ys > s 


C(t,s\r) = 


Hi 


F r (t + s) + 


2 Tr 


\A 


dr g(r)F r (t + s — 2r) 


2 T r 


Hi 


2 Tr 


d t 5{r)F r {t + s — 2r) 


c JO 


s—yoo 


y/g{t)g(s) 

H\F r (t + s) 

\fg(t)g{s) \fg{t)g{s) 

+2 T c (ts) {2 ~ d)/A [ drr d/2 ~ 1 F r (t + s- 2r) 

Jo 

s i-i/2 2 T c y(2-j)/4 f d'u u d ! 2 ~ l (47T [(y + 1) — 2'«])~ rf/2 exp 


(A.6) 


S 7/ + 1 — 2u 


where in the second line the singular contribution in g(r) from (IA.3[) has been explicitly intro¬ 
duced to demonstrate the cancellation of the LA-dependent terms. In the third line, the regular 
large-time asymptotics of g(t ) has been inserted and finally, F r (t) is expanded for large times. 
For d > 2, we now have g(t ) ~ f° instead of t d F~ 1 and an analogous calculation gives 


C(t,s-,r) = s 1 d / 2 2T c / d uu d ^ 2 1 (47t [(y + 1) — 2u}) d ^ 2 


exp 


s y + 1 — 2u 


(A.7) 


which proves the assertion in section 4. For r = 0, the two-time autocorrelator and its scaling 
function (14.10(1 is readily obtained. 

Finally, we examine the special case d — 2. Indeed, the treatment is analogous to the 
one of the AD magnetic spherical model J74] [79]. From (lA.lj) . one readily Ends g(p). For 
large times, it can be shown that g(t) — —Hi/(2T c )S(t) + < 7 2 /lnf [SI IS, 109] . The spherical 
parameter g(t) enters into the two-time response and correlation functions only through ratios 
g(t)/g(s ) ~ Inf/Ins = ln(ys)/ In s ~ 1 + In y/ Ins. Therefore, these logarithmic terms only 
arise as additive logarithmic corrections to scaling, but do not affect the scaling form itself 
[74] [79]. This can also be seen in the explicit results (14.71) . (14.101) . since the scaling functions are 
continuous in d, even at d — 2. This is different, however, for the interface width, where the 
amplitudes eq. (I4.11|) diverge for d —* 2. This case must therefore be analysed separately. In 
analogy with the treatment which led to (IA.4D for d < 2, we have first the cancellation of the 
LA-dependent terms and then estimate the large-f behaviour of the remaining integral 

w 2 {t) ~ 2T C In f J dr ^—F 0 (2t - 2r) = 2T c lnf C~ l (h/ 7 )) ^ (^o(2r)) (p)^ (f) 

~ -^-4- hi t (f) ~ Inf (A.8) 


Again, the non-universal amplitude g 2 cancels. We recognise the integral as a Laplace con¬ 
volution and estimate the leading small-p behaviour of the two factors. First, we recall 
(£(l/lnr)) (p) rs./ — \p(Ce + lnp)] 1 , see [74] [79]. The second factor is the Laplace trans¬ 
form of F 0 (2t) in 2D , which is evaluated in analogy with similar integrals treated above, as 
follows (with a cut-off p; E\ is the exponential integral [HSi eq. (5.1.11)]) 


(«(2r)) (p) 



dr e pr (87rr) 1 


87T 


Eiipn) 


p —>o 


C E + lnp 
8n 
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As in ca, terms only depending on In q must be absorbed into the ‘finite’ contribution coming 
from the integral JJ 7 , which leads here to non-leading terms. This proves the middle line in 
gTU). 


A.2 Ageing at low temperatures 


For T < T c , define the auxiliary variable m 2 := 1 — T/T c . First, we find g(p ) from (14.11) . with 
f(p) explicitly given in (IA. II) . This gives 


9 ip) - 


iA J_ 

2 T c m? 

iAJ_ 

2 T c m? 


_l_ H\A d d/2 I 
+ dm* P + • • • 

+ IA(o)p + + • • • 


if 0 < d < 2 
if 2 < d < 4 


(A.9) 


where we used T c + T/m 2 = T c /m 2 . Further regular terms appear for d > 4, 6,..., but will only 
give rise to corrections to the leading scaling behaviour. This gives in turn 


g(t) 


StsM ( t ) + ;&( 8,r )~“ /2 *~ 1 ~ i/2 + '-- ; ifo < d < 2 

+ + + ; if 2 < d < 4 


(A.10) 


The response function is once more obtained by straightforward asymptotic expansion, 
which leads again to (j4.6|) with the scaling function fn{y) being now given by (14.181) . 

Next, we analyse the interface width. The techniques to be used are quite close to the one 
applied in the critical case, but the leading term turns out to be of a different form (let d < 2 
for simplicity) 


w 2 (t) = 


H i F o(2t) +2T f dr 4^F 0 (2t-2r) 


9(t) 


9(t) 




i t \ r 1 - / r \-i-d /2 

git) \ 1 + ^T c ) Fo{2t) + 2T J 0 dT (?) F °( 2t ~ 2T ) 

~ 4 m 2 t + 2Tt l+d / 2 C ~ 1 ((£(r- 1 - d/2 )(p)) (£F 0 (2r)(p))) (t) 

= 4m 2 f + 2Tt 1+d/2 1 ( d/ 2 ~ d/2) (^~ 1 /~ 1 ) (t) 

= 4m 2 t + O (Tt 1_d/2 ) 


(A.ll) 


In the second line, we used both the singular and the leading regular term for g{r). In contrast 
to the critical case, the singular contribution does not cancel with the other H \-dependent 
term, but combines into a new, leading term which depend on the ratio T/T c and vanishes 
when T -A- T c . Therefore, the leading long-time behaviour only depends on the ratio T/T c 
and the short-ranged nature of the initial correlations. The ‘thermal’ convolution term can be 
calculated via a Laplace convolution as before, but turns out to provide merely a correction to 
scaling. For d > 2, the leading correction is O (T) and the further derivatives of the 5-function 
only generate, via partial integrations, further sub-leading corrections to scaling. Hence the 
leading result w 2 {t ) ~ 4m 2 t, valid for t large enough, is independent of d, as asserted in (14.21)1 . 

The two-time autocorrelator is analysed in the same way (assume d <2 and let y = t/s > 1) 
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C(M) = 


H 1 F 0 (t + s) 
y/9(t)g(s) 


2T / dr 


SCO 


m 4 r(— d/2)d (ts ) ( ' 2+d ^ 4 

A ^ tt )'*/ 2 (t + s ) d / 2 


1 + 


\/9{t)g{s) 
1 T 


m 2 T r 


+ 


F 0 (t + s — 2r) 
2 T 


(47r) d / 2 


= 2 d+2 m 2 s 7/ (2+ci)/4 (7/ + l) _d/2 + O (Ts 1_<i/2 ) 


( ts )( 2 +( i )/ 4 / dr r - i - d / 2 ( t + s _ 2r )- d / 2 

Jo 

(A.12) 


as stated in (j4.19|) . For d > 2, the leading correction is O(T). 

Since all results are continuous in d at d — 2, a separate analysis of the 2.D case is not 
necessary. 


Appendix B. On the Yeung-Rao-Desai inequalities 

The analogue of the well-known YRD-inequality [104] , originally formulated for ageing magnetic 
systems, is: for a growing interface with a non-conserved dynamics and an uncorrelated initial 
state, the autocorrelation exponent Xq satisfies the bound 

A c>^(d + zb) (B.l) 

Such bounds may serve as checks on exponents, estimated from simulational or experimental 
data (i.e. for phase-ordering magnets with T < T c , one has 6 = 0, hence Ac > d/2 [ 104] ). 

To see this, recall first the argument [ 104] to bound the two-time height autocorrelator 

C(t, s)= Idk C(t , s; k , -k) < Jdk \Jd(t, k ) d{s, k) (B.2) 

Eq. (1.2) gives the scaling form of the single-time correlator C(t, k ) = C(t,t ; k, —k ) 

C(t, k) = -^= [dr e~ ikr t~ b F c (l; rF 1 ^) = L{t) d ~ zb C ( kL(t )) 

W Jv 

where the typical length scale L(t) ~ t x F and C is a scaling function. In order to estimate the in¬ 
tegral, recall that the height fluctuations 5h(t, k ) and 5h(s, —k ) in C(t , s; k , —k) should become 
uncorrelated over distances larger than A L > 27ra/|fc| such that C(t,s ; k , — k) —y 0 rapidly if 
(L(t)-L(s))\k\ > i m- For uncorrelated initial conditions, one has lim s ^o C*(s, fc) ~ 
limfe-^-o C(s, k ) = 0(1). In the scaling limit with y = t/s large enough, the bound (IB.2|) gives 

lim C(t, s) ~ L(t)~ Xc 

t,s—too 

/» 27 ra / L(t) /* 27 ra 

< L{tY d ~ zb)/2 C d / dkk^C 1 ' 2 (kL{t)) = L{t)- [d+zb)/2 C d / duu i_1 c 1/2 (ti) 

Jo Jo 

{C d is a constant) and (IB.ID follows. 

| | _ 

Alternatively, if the waiting time s is itself already in the scaling regime, one has C(s , k) ~ 
fczh-d w j 1 i c j 1 leads to Ac > zb, with an magnetic analogue already given in [104] , 
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